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distribution of an unordered eigenvalue and certain random determinant properties. Based on the analytical 
capacity expressions, we investigate the impact of the system and channel characteristics, such as the 
antenna configuration and the relay power gain. We also demonstrate a number of interesting relationships 
between the dual-hop AF MIMO relay channel and conventional point-to-point MIMO channels in various 
asymptotic regimes. 
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I. Introduction 

The relay channel, first introduced in [1,2], has been considered in recent years as a means to improve 
the coverage and rehabihty, and to reduce the interference in wireless networks [3-1 1]. Generally speaking, 
there are three main types of relaying protocols: decode-and-forward (DF), compress-and-forward (CF), 
and amphfy-and-forward (AF). Of these protocols, the AF approach is the simplest scheme, in which case 
the sources transmit messages to the relays, which then simply scale their received signals according to a 
power constraint and forward the scaled signals onto the destinations. 

Point-to-point multiple-input multiple-output (MIMO) communication systems have also been receiving 
considerable attention in the last decade due to their potential for providing Unear capacity growth and 
significant performance improvements over conventional single-input single-output (SISO) systems [12, 
13]. Recently, the application of MIMO techniques in conjunction with relaying protocols has become a 
topic of increasing interest as a means of achieving further performance improvements in wireless networks 
[14-18] 

In this paper we investigate the ergodic capacity of AF MIMO dual-hop systems. This problem has 
been recently considered in various settings. In [19], the ergodic capacity of AF MIMO dual-hop systems 
was examined for a large numbers of relay antennas K, and was shown to scale with log K. Asymptotic 
ergodic capacity results were also obtained in [20] by means of the replica method from statistical physics. 
In [21, 22], the asymptotic network capacity was examined as the number of source/desination antennas M 
and relay anteimas K grew large with a fixed-ratio K/M (5 using tools from large-dimensional random 
matrix theory. It was demonstrated that for /3 — >^ oo, the relay network behaved equivalently to a point- 
to-point MIMO link. The results of [21,22] were further elaborated in [23] where a general asymptotic 
ergodic capacity formula was presented for multi-level AF relay networks. Recently, the asymptotic mean 
and variance of the mutual information in correlated Rayleigh fading was studied in [24]. All of these 
prior capacity results, however, were derived by employing asymptotic methods (i.e. by letting the system 
dimensions grow to infinity). To the best of our knowledge, there appear to be no analytical ergodic 
capacity results which apply for AF MIMO dual hop systems with arbitrary finite anteima and relaying 
configurations. 

In this paper we derive new exact analytical results, simple closed-form high SNR expressions, and 
tight closed-form upper and lower bounds on the ergodic capacity of AF MIMO dual-hop systems. In 
contrast to previous results, our expressions apply for any finite number of MIMO anteimas and for 
arbitrary numbers of relay antennas. The results are based heavily on the theory of finite-dimensional 

random matrices. In particular, our exact ergodic capacity results are based on a new exact expression 
which we derive for the exact unordered eigenvalue distribution of a certain product of finite-dimensional 
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random matrices, corresponding to the equivalent cascaded AF MIMO relay channel. In prior work [22], an 
asymptotic expression was obtained for this unordered eigenvalue density. However, that asymptotic result, 
which serves as an approximation for finite-dimensional systems, was rather complicated and required the 
numerical computation of a certain fixed-point equation. Our result, in contrast, is a simple exact closed- 
form expression, involving only standard functions which can be easily and efficiently evaluated. In addition 
to the unordered eigenvalue distribution, we also present a number of new random determinant properties 
(such as the expected characteristic polynomial) of the equivalent cascaded AF MIMO relay channel. These 
results are subsequently employed to derive simplified closed-form expressions for the ergodic capacity in 
the high SNR regime, as well as tight upper and lower bounds. Again, these random determinant properties 
are exact closed-form analytical results which apply for arbitrary antenna and relaying configurations, 
and are expressed in terms of standard functions which are easy to compute. As a by-product of these 
derivations, we also present some new unified expressions for the expected characteristic polynomial and 
expected log-determinant of semi-correlated Wishart and pseudo-Wishart random matrices. 

Based on our analytical expressions, we investigate the effect of the different system and channel param- 
eters on the ergodic capacity. For example, we show that when either the number of source, destination, 
or relay antennas, or the the relay gain grows large, the AF MIMO dual-hop capacity admits a simple 
interpretation in terms of the ergodic capacity of conventional single -hop single-user MIMO channels. In 
the high SNR regime, we present simple closed-form expressions for the key performance parameters — the 
high SNR slope and the high SNR power offset — which reveal the intuitive result that the multiplexing gain 
is determined by the minimum of the number of antennas at the source, destination, and relay, whereas the 
power offset is a more intricate function which depends on all three. For example, we show that by adding 
more antennas at the destination, whilst keeping the number of source and destination antennas fixed, may 
lead to a significant improvement in the high SNR power offset; however the relative gain becomes less 
significant as the initial number of destination antennas is increased. Our analytical expressions also reveal 
the interesting result that the ergodic capacity of AF MIMO dual-hop channels is upper bounded by the 
capacity of a SISO additive white Gaussian noise (AWGN) channel. 

The remainder of this paper is organized as follows. Section |II] presents the AF MIMO dual-hop system 
model under consideration. Section JII] presents our new random matrix theory contributions, which are 
subsequently used to derive the exact, high SNR, and upper and lower bound expressions for the ergodic 
capacity in Sections |IV] and jV] Section |Vl] summarizes the main results of the paper. All of the main 
mathematical proofs have been placed in the Appendices. 



3 




Source Relay Destination 



Fig. 1. Schematic diagram of a MIMO dual-hop system, where there is no direct Unk between source 
and destination. 

II. System Model 

We employ the same AF MIMO dual-hop system model as in [21,22]. In particular, suppose that there 
are source antennas, rir relay antennas and destination antennas, which we represent by the 3-tuple 
{ns,nr,na). All terminals operate in half-duplex mode, and as such communication occurs from source 
to relay and from relay to destination in two separate time slots. It is assumed that there is no direct 
communication link between the source and destination, as sketched in Fig. \T\ The end-to-end input-output 
relation of this channel is then given by 

y = H2FH1S + H2Fn„^ + (1) 

where s is the transmit symbol vector, and are the relay and destination noise vectors respectively, 
F = \/a/ {ur (1 + /o))Ira,. (a corresponds to the overall power gain of the relay terminal) is the forwarding 
matrix at the relay terminal which simply forwards scaled versions of its received signals, and Hi G (j^^x^s 
and H2 G c^dxn^ denote the channel matrices of the first hop and the second hop respectively, where their 
entries are assumed to be zero mean circular symmetric complex Gaussian (ZMCSCG) random variables of 
unit variance. The input symbols are chosen to be independent and identically distributed (i.i.d.) ZMCSCGs 
and the per antenna power is assumed to be p/n^, i.e., E {ss^"} = {p/ng) In^- The additive noise at the 
relay and destination are assumed to be white in both space and time and are modeled as ZMCSCG with 
unit variance, i.e., E |n„^nn^| = and E |n„^nJi^| = We assume that the source and relay have 
no channel state information (CSI), and that the destination has perfect knowledge of both H2 and H2H1. 

The ergodic capacity (in b/s/Hz) of the AF MIMO dual-hop system described above can be written as 
[20-22] 

C7 = ^i?{log2det(l + R,R-i)} (2) 



where and R„ are x matrices given by 



and 



Rs — — H2HiH|h1 



(3) 



Rn. — + aH2H2 



(4) 



respectively, with 



Using the identity 



a 



n-r (1 + p) 



(5) 



det (I + AB) = det (I + BA) , 
Q can be alternatively expressed as follows 

C (p) = |log2 det (in^ + ^h1h|r-1H2Hi^ | . 
Next, we utilize the singular value decomposition to write H2 = U2D2V2, where 

D2 = diag { Ai , . . . , Ainin(nd,n,) } 



(6) 



(7) 



(8) 



is an nd x diagonal matrix, with diagonal elements pertaining to the increasing ordered singular values, 
and U2 G C'^dXna g Qn^xn-r ^j-g unitary matrices containing the respective eigenvectors. Since Hi 

is invariant under left and right unitary transformation, the ergodic capacity in dTJl can be further simplified 



as 



C ip) = \e |log2 det (in^ + ^H^^Hi 
2 V n.o 



where 



* = < 



diag<' ^^,...,^^;^,0^_.^[> , nr>nd. 



It is then easily established that 



C (p) = ]-E |log2 det f + ^hIlHi 
2 V n,. 



(9) 



(10) 



(11) 



where h| ~ CJ\fn,,q (0, Ig), with q = min (n^, n^), and 

L = diag{A2/(l + aA?)}^^^ . (12) 

Equivalently, we can now write 

C (p) = ^ ^°S2 (l + a (A) d\ (13) 

where s = min(ns,g), A denotes an unordered eigenvalue of the random matrix h|LHi, and /a (•) 
denotes the corresponding probability density function (p.d.f.)- Although the distribution of A has been 
well-studied in the asymptotic antenna regime [21, 22], currently there are no exact closed-form expressions 
for /a(-) which apply for arbitrary finite-antenna systems. 

III. New Random Matrix Theory Results 

In this section, we derive a new exact closed-form expression for the unordered eigenvalue distribution 
/a(-) of the random matrix h|lHi. We also present a number of other key results, such as random 
determinant properties, which will prove useful in subsequent derivations. It is convenient to define the 
following notation: = A^, /3j = A|/ (l + aXf) {i = 1, . . . ,q), and p = max {n^, rir). 

To derive the unordered eigenvalue distribution /a(-)> we first need to establish some key preliminary 
results, as given below. 

Lemma 1: The marginal p.d.f. of an unordered eigenvalue A of h|LHi, conditioned on L, is given by 

1 J-^ )^ns+k-q-l^-X/f3igq-ns-l 

where Di f^ is the A;)th cofactor of a. q x q matrix D whose (m, n)th entry is 

{D}^_„ = (15) 
Proof: See Appendix II-AI □ 
This lemma presents a new expression for the unordered eigenvalue distribution of a complex semi- 
correlated central Wishart matrix. In prior work [25], two separate alternative expressions for this p.d.f. 
were obtained for the specific scenarios Us < q and Ug > q respectively; the latter casqj being a complicated 
expression in terms of determinants with entries depending on the inverse of a certain Vandermonde matrix. 
Here, Lemma [7] presents a simpler and more computationally-efficient unified expression, which applies 
for arbitrary and q. 

To remove the conditioning on L in Lemma [H it is necessary to establish a closed-form expression for 

'For this case (n^ > g), the random matrix H^LHi has reduced rank and the corresponding distribution, conditioned on L, is 
commonly referred to as pi-etfiio- Wishart [26]. 



the joint p.d.f. of • • • , /3g. We will also require the p.d.f. of an arbitrarily selected /? G • • • ,Pq}. 
These results are given in the following lemma. 

Lemma 2: The joint p.d.f. of {0 < /3i < • ■ • < /3g < 1/a} is given by 



1 pf~ie~A 



where 



/C 



{J^.^^'^{q-i + '^)'^{p-i + ^)) ' ■ (IV) 



The p.d.f. of an unordered (randomly-selected) /? S • • • ,/3g} is given by 
where 

Proo/- See Appendix II-BI □ 
Having established the results in Lemma [7] and Lemma \2l we are now ready to derive the desired 
unconditional unordered eigenvalue distribution /a(-)' given below. 

Theorem 1: The marginal p.d.f. of an unordered eigenvalue A of h|lHi is given by 

h (A) = ^ E E E T\n a + k) A^^-^"^^-^-^^/^ W^-i (2^) G,. (20) 

1=1 k=q-s+l i=0 V s y -T ; 

where (•) is the modified Bessel function of the second kind and Gi^k is the (Z, A;)th cofactor of a g x g 
matrix G whose (m, n)th entry is 

{G}„,„ = a^-P-^-^+^r (p-g + m + n-l);7(p-g + m + n-l,p + g, 1/a) (21) 

with [/ (•, •, •) denoting the confluent hypergeometric function of the second kind [27, Eq. 9.211.4]. 

Proof: See Appendix II-CI □ 
We note that an asymptotic expression for /a(-) has been considered previously in [22], based on large- 
dimensional random matrix theory. However, that asymptotic p.d.f. result, which serves as an approximation 
for finite-dimensional systems, is not in closed-form, requiring the numerical computation of a certain 
fixed-point equation. Indeed, to further facilitate computation of the asymptotic eigenvalue p.d.f. in [22], 
an algorithmic approach with certain heuristic elements was also presented. Our result in Theorem [7] in 
contrast, gives the exact eigenvalue p.d.f. which applies for arbitrary finite system dimensions, and is 
presented in a simple closed-form involving only standard functions which can be easily and efficiently 
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evaluated. In the following section, this result will be employed to evaluate the ergodic capacity of AF 
MIMO dual-hop channels. 

Corollary 1: For the special case (1, 1, 1), the unordered eigenvalue p.d.f. ( [20l ) reduces to 



fx 'W = 2e 



VXKi (iVx) + Ko (2Vx) 



(22) 

Proof: The proof is straightforward and is omitted. □ 
We note that this special case has also been derived previously in [28]. 

Corollary 2: Let L = diag {A^}J_j^. Then, the marginal p.d.f. of an unordered eigenvalue A of H{LHi 
is given by 

or * \(n.+2fc+p+«-2g-3)/2 

hW^-T. E r(n.-, + ;c) W--(2v/A)g,. (23) 

1=1 k=q-s+l ^ < / 

where G/ ^ is the {I, k)th cofactor of a q x q matrix G whose (m, n)th entry is 

{G}^^^ = T{p-q + m + n-l). (24) 
Proof: The result is obtained by taking the limit as a ^ in (l20l) . □ 
This result will be used to study the capacity of AF MIMO dual-hop channels in the high SNR regime. It 
is also worth noting that (1231 ) can be applied to the ergodic capacity analysis of Rayleigh-product MIMO 
channels [29,30]. 

Fig. |2] compares the analytical result presented in Theorem |7] with Monte Carlo simulations. We plot 
the p.d.f. of the unordered eigenvalue A with system configuration (2,3,4). The simulated p.d.f. curve is 
based on 100,000 channel realizations. The figure shows that the analytical result is in agreement with the 
simulations. 

Fig. [3] shows the analytical result presented in Theorem [7] and Corollary \2\ The curves corresponding 
to p = dB, p = 10 dB, and p = 20 dB are generated using ( [23] ) while the "Rayleigh Product" curve 
is generated using ( |23l ). We can see that the exact unordered eigenvalue distribution converges to the 
unordered eigenvalue distribution of the Rayleigh product channel as a — > cxd, as expected. 

Fig. |4] compares our exact unordered eigenvalue distribution, based on ( |20l ). with the corresponding 
asymptotic eigenvalue distribution presented in [22], for the random matrix HjLHi/(nsnr) with different 
system configurations. We use the same simulation parameters as in [22, Fig. 5 (a)], setting a = l/rir 
and Ur/us = 1/2. We clearly see the convergence of the exact and asymptotic p.d.f.s as the numbers of 
antennas become large (eg. the (16, 8, 16) scenario), however when the systems dimensions are not so large 
(eg. the (2,1,2) and (4,2,4) scenarios), the asymptotic eigenvalue p.d.f. exhibits noticeable inaccuracies 
with respect to our new exact result in (l20l) . 

The following theorems present new closed-form random determinant properties, involving the random 
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Fig. 2. Comparison of the analytical and Monte Carlo-simulated unordered eigenvalue p.d.f. of H{LHi. 
Results are shown for (2, 3, 4) system configuration, with a = 2. 
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Fig. 3. Comparison of the analytical unordered eigenvalue p.d.f. of h|lHi and h|lHi for different p. 
Results are shown for a (2, 3, 4) system configuration, with a = 2. 
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Fig. 4. Comparison of the analytical unordered eigenvalue p.d.f. of H|LHi/(nsnr) for different system 
configurations. Results are shown for a = l/rir and rir/us = 1/2. 



matrix h|lHi. These results will be applied to derive tight bounds on the ergodic capacity. 
Lemma 3: The expected determinant of + {pa/ng) H{LHi, conditioned on L, is given by 

det (A) 



£; <i det ( + — hJlH] 



where A is a g x q matrix with entrie: 



I 



(25) 



! n<q- Us, 

/5m"' fl + ^/?mK-(? + n)), n>q-ns. 



(26) 

□ 



Proof: See Appendix II-DI 
This theorem presents a new expression for the expected characteristic polynomial of a complex semi- 
correlated central Wishart matrix. In prior work [31, 32], alternative expressions were obtained via a differ- 
ent approach (i.e. by exploiting a classical characteristic polynomial expansion for the determinant). Those 
results, however, involved summations over subsets of numbers, with each term involving determinants of 
partitioned matrices. In contrast, our result in Lemma |7] is more computationally-efficient, involving only a 
single determinant with simple entries. Moreover, it is more amenable to the further analysis in this paper, 
leading to the following important theorem. 



'When 



q < Us, {A}^,„ = /3r ' (l + {ns-q + n)) . 
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Theorem 2: The unconditional expected determinant of + (pa/ng) h|lHi is given by 

E |det + ^hJlHi j I = /C det (S) (27) 



where H is a g x g matrix with entries 

171,11 

with T=p — q + m + n, and 



'^i?T-_i(a), n < q — Us 



(28) 



7?,(a) =r(r)[/(r,^j + g,l/a) . (29) 
Proof: Utilizing Lemma |21 [33, Lemma 2] and (1112b yields the desired result. □ 
Lemma 4: Let 



* = s : _ , (30) 

[ LHiHj, q<ns. 
The expected log-determinant of conditioned on L, is given by 

i det(Yfc) 

E {In det (*)| L} = ^ (n. - s + fe) + ^^^^^^^^^ -— (31) 

where (•) is the digamma function [27], and Yjt is a g x g matrix with entries 

{Y.},..,,. = -i (32) 

/5;;j Mn/3m, n = k. 



When q = s, (1311) reduces to 

s 

S{lndet(*)|L} = ^VK-s + fc) + lndet(L) . (33) 

k=l 

Proof: See Appendix II-EI □ 
We note that the above expected natural logarithm of the determinant for q > Ug has been investigated 
in [34], where the derived expression is rather complicated, involving summations of determinants whose 
elements are in terms of the inverse of a certain Vandermonde matrix. We also note the q < Ug and 
q = Ug = s cases have been considered in [32, 35]. Our result, in contrast, gives a simple unified expression 
which embodies all of these cases. Moreover, based on Lemma |?1 we obtain the following important 
theorem. 

Theorem 3: The unconditional expected log-determinant of $ is given by 

s q 

^{lndet(*)} = ^V'K -s + A;) +/C ^ det(Wfe) (34) 

fc=l k=q—s+l 
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where Wjt is a q x q matrix with entries 



{W.}_ ={ ^ (35) 

Wn(a)) n = k 



where r and are defined as in ( |29l ). and 

i=0 V ^ / \ ^""^ 



where gi{-) denotes the auxiliary function 

<7;(x) = e"^i+i(x) (37) 

with (•) denoting the exponential integral function of order I + 1. 
When q = s, (l34l) reduces to 

q-l i 2j 2g-l-2 z.. _ / _ 

i5;{lndet(*)} = j;^(n,-s + A:) + j;j;^ E ( A; M J' 9) 

k=l i=0 j=0 1=0 k=0 ^ ^ 



^ m=0 ^ 



Proof: See Appendix II-FI □ 

IV. Ergodic Capacity Analysis 

In this section we present new analytical expressions for the ergodic capacity of AF MIMO dual-hop 
systems. 

A. Exact Expression for Ergodic Capacity 
Substituting into ([T3]) we obtain 

q q q+ris-l /q+n,~l\^q+n,~l-i 

cM^icY, T. E V (39) 

1=1 k=q-s+l i=0 \ s 'i • ; 



where 



(40) 



The integral in (l40b can be evaluated either numerically, or can be expressed as an infinite series involving 
Meijer-G functions. These results are confirmed in Fig. |5l where we compare the exact analytical capacity 
of AF MIMO dual-hop systems, based on (l39l) and (l40l) . with Monte-Carlo simulated curves for two 
different antenna and relay configurations. In both cases, there is exact agreement between the analysis 
and simulations, as expected. 
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1) Analogies with Single-Hop MIMO Ergodic Capacity: Let C^^~^™^ {ns,n^, p) denote the ergodic 
capacity of a conventional single -hop i.i.d. Rayleigh fading MIMO channel matrix H G Qn^xn^^ ^^^^ 
transmit and receive antennas, and average SNR p; i.e. 

^SH-MiMO(„^^ p^^E |log2 det (^I„,, + -^HHt) I . (41) 

Here, we demonstrate four particular cases for which the AF MIMO dual-hop channel relates directly to 
single-hop MIMO channels, in terms of ergodic capacity. 

• As the number of relay antennas grows large, i.e. — > oo, the ergodic capacity of AF MIMO 
dual-hop systems becomes 

lim C ip) = Ic^H-MiMO L n,, . (42) 

A proof is presented in Appendix III-AI Note that a similar phenomenon has been derived in [19], 
for the special case Ug = na. Here, (l42l) generalizes that result for arbitrary source and destination 
antenna configurations. 

• As the number of source antennas grows large, i.e. oo, the ergodic capacity of AF MIMO 
dual-hop systems becomes 

lim C ip) = IC^H-MIMO (^^^ „^ _ 1 ^SH-MIMO f^^^^^ »\ (43) 
n^^oo 2 2 \ 1 + p J 

A proof is presented in Appendix III-BI Interestingly, we see that as p grows large, the right-most 
term in ( |43] ) disappears, and the AF MIMO dual-hop capacity becomes equivalent to one half of the 
ergodic capacity of a single -hop MIMO channel with rir transmit antennas, receive antennas, and 
average SNR a. 

• As the number of destination antennas grows large, i.e. rid 00, the ergodic capacity of AF MIMO 
dual-hop systems becomes 

lim C7(p) = ^CSH-MiMO(^^^„^^^^) (44) 

rid— ►oo 2 

The result is trivially obtained by directly taking — > oo in (fTTI) . We see that the AF MIMO dual-hop 
capacity becomes equivalent to one half of the ergodic capacity of a single-hop MIMO channel with 
Us transmit antennas, rir receive antennas, and average SNR p. 

• As the power gain of the relay grows large, i.e. a ^ 00, the ergodic capacity of AF MIMO dual-hop 
systems becomes 

lim C(p) = ic7^^^-^i^OK,g,p) . (45) 

a— >oo 2 
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The result is trivially obtained by directly taking a — > oo in ([TT] ). Thus we see the interesting result 
that even as the relay power gain becomes very large, the capacity of AF MIMO dual-hop channels 
remains bounded, and in fact becomes equivalent to one half of the ergodic capacity of a single-hop 
MIMO channel with Ug transmit antennas, q = min(nr,n(i) receive antennas, and average SNR p. 
We note that for each of the cases (l42l)-(|45]). closed-form expressions can be obtained by directly 
invoking known results from the single-hop MIMO capacity literature (eg. see [31]). 

In order to obtain further simplified closed-form results, it is useful to investigate the ergodic capacity 
in the high SNR regime. This is presented in the subsection below. 

B. High SNR Capacity Analysis 

For the high SNR regime, we consider two important scenarios; namely, one where the source and relay 
powers grow large proportionately, and one where the source power grows large but the relay power is 
kept fixed. 

1) Large Source Power, Large Relay Power: Here we have a ^ oo, p ^ oo, with a/ p = (3, for some 
fixed j3. Then pa ^ ^ and a P/ur, and the ergodic capacity at high SNR reduces to 

C {p)\a,p^oo,a/p=p = \e {log2 det + } (46) 

where L = diag {Af / (l + (/3/n^) Af) We can express (l46l ) in the general form [34] 

C (p)L,p^oo„a/p=/3 = ^oo - £00) + o (1) (47) 

where 3dB = 101ogiQ(2). Here, the two key parameters are S^o, which denotes the high-SNR slope in 
bits/s/Hz/(3 dB) given by 

S^= lim ^^^^'"■P-°o."/P=/^ (48) 

Q,p-^oo log2(p) 

and Coo, which represents the high-SNR power offset in 3 dB units given by 

£00 = lim log2(/o) 5 . (49) 

a,p^oo y boo J 

From (l46l) . we can evaluate Soo and £00 in closed-form as follows. 

Theorem 4: For the case a ^ 00, p ^ 00, with aj p = /?, the high-SNR slope and high-SNR power 
offset of AF MIMO dual-hop systems are given by 

5oo = | bit/s/Hz/(3dB) (50) 
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anil 



£oo{ns,nr,nd) = log 



sin 2 



^2^^ (r^s + k - s) + IC det(Wfc) 

k=l k=q—s+l 



(51) 



respectively, where is a. q x q matrix with entries 

I- "'J m,n 



(52) 



n = k. 



For the case q = s (i.e. corresponding to min{ns,nr,nd) = rid or min{ns, nr,nd) = rir), the high SNR 
power offset (ISTl ) admits the alternative form 



Coo{ns,nr,nd) = log 



nsUr 



1 



p 



sln2 

2q-l-2-k 



2q-l-2 

k 



E^in.-s+k) + j:EE E 

k=l 4=0 j=0 1=0 k=0 

(p+q-k-2 
^{p + q-k-1)- 9 
m=0 



(53) 

□ 



Proof: See Appendix III-CI 
Interestingly, we see that the high SNR slope depends only on the minimum system dimension, i.e. s = 
min(ns, n^, nrf), whereas the high SNR power offset is a much more intricate function of n^, n^, and n^. 
Fig. |5] depicts the analytical high SNR capacity approximations for AF MIMO dual-hop systems, based 
on dSOt and (ISTT ). These approximations are seen to converge to their respective exact capacity curves for 
quite moderate SNR levels (eg. < 20dB). 

It is important to note that Theorem |?] presents an exact characterization of the key high SNR ergodic 
capacity parameters, S^o and £oo(')' arbitrary numbers of antennas at the source, relay, and destination 
terminals. We now examine some particularizations of Theorem |?1 in which these expressions reduce to 
simple forms. 

Corollary 3: Let Ur = 1. Then Soo = 1/2, and £oo(') reduces to 





1 


(7) 




= Inn^ 


+ o{l) 



(54) 



ns — > 00, and as such we have 



lim_ Cooiris, 1, rid) = log2 ( ^) ~ f"^ 



(55) 



''Note that here we explicitly indicate the dependence of the high SNR power offset on n^, n^, and n^. 
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TABLE I. High SNR offset as function of rid, where = 2, = 3 and (3 = 2 



rid 


4 


6 


8 


10 


12 


14 


Coo (dB) 


2.593 


1.573 


1.147 


0.88 


0.73 


0.622 


offset as : 


unction of rir, where = 2, = 4 anc 


1/3 = 2 




3 


5 


7 


9 


11 


13 


Coc (dB) 


2.593 


1.251 


0.847 


0.636 


0.493 


0.429 



Corollary 4: Let Ud = 1. Then Soo = 1/2, and £oo(') reduces to 



loe 



/3 



1 



m=Q 



(56) 



m=0 VP/ 



(57) 



le relative power gain factor /3 on the ergodic 



In this case, as Us grows large we have 

lim £oo(ns,nr, 1) = log2 ( ^ ) - 7^ 
Based on these results, we can easily examine the effect of t 

capacity. In particular, noting that gi (x) in (|37] ) is a monotonically decreasing function of x in the interval^ 
[0,oo), we see that increasing (3, whilst having no effect on the high SNR capacity slope Soo, results in 
decreasing the high SNR power offset £oo(-)' ^i^d therefore increasing the ergodic capacity in the high 
SNR regime. 

Corollary 5: Let Ug = rir = 1- Adding k destination antennas, while not altering 5oo, would reduce the 
high SNR power offset as 

S{nd, k) = Coo (1, 1, rid + k) - Coo (1, 1, rid) 

l=nd 

Note that, to obtain this result, we have invoked the definition of the digamma function [27]. Since 
Qi {x) > for X e [0,cxd), it is clear that the high SNR power offset £oo(') in ( l58l ) is a decreasing 
function of k, thereby confirming the intuitive notion that adding more antennas to the destination terminal 
has the effect of improving the ergodic capacity. 
Example 1: With respect to 13 = 1, 



£oo (1,1, 2) = £00 (1,1,1)- 2.58 dB (59) 
£00 (1,1, 3) = £00 (1,1,1)- 3.46 dB (60) 
£00(1,1,00) = £00(1,1,1) -5.08 dB (61) 

where £00 (1,1,1) = 7.57 dB. 

"'This conclusion is easily established by noting that A/Ax {gi (x)) = [Ei+i (x) — Ei (a;)], and using [36, Eq. 5.1.17]. 
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Up 
12 - 




Fig. 5. Comparison of exact analytical, high SNR analytical, and Monte Carlo simulation results for 
ergodic capacity of AF MIMO dual-hop systems with different antenna configurations. Results are shown 
for a/p = 2. 




Fig. 6. High SNR power offset shift, in decibels, obtaining by adding either (a) one antenna to the 
destination, (b) two antennas to the destination, or (c) four antennas to the destination. Results are shown 

for Us = Ur = I and a/p = 2. 
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Fig. [6] illustrates the relationship in Corollary\5\ where the high SNR power offset shift 5{n(i, k) is plotted 
against n^, for k = I, k = 2, and A; = 4. As expected, for a fixed value of k, 6{n(i,k) is an increasing 
function of na, approaching a limit of dB as oo. Table I and Table II present the high SNR power 

offset as a function of and respectively, for Ug = 2. We see that when (resp. n^) is small, then 
a small increase in n^; (resp. rir) yields a significant improvement in terms of the high SNR power offset. 
However, in agreement with Fig. |6l adding more and more antennas yields diminishing returns. 

2 ) Large Source Power, Fixed Relay Power: Here we take /? — > oo and keep a fixed. Then, noting that 
po-lp^oo ~* ct/rir, the ergodic capacity reduces to 

lim C ip) = |log2 f 1 + I (62) 

where A denotes an unordered eigenvalue of Hj^LHi. Using Corollary^ we can evaluate this constant as 

limC(p) = — V V Gi,k:Fi,k (63) 

1=1 k=q-s+l 



where 



Ti,, = r In (l + ^y] y("=+2/c+p+^-2,-3)/2^^^^_^^_^ ^^^^ (64) 

JO V ngUr / 



To evaluate the remaining integral in (I621 ). we first express the logarithm in terms of the Meijer G-function 
as [37, Eq. 8.4.6.5] 



( " A 


1, 








:) 


\ risHr 


1, 





(65) 



and then apply the integral relationships [27, Eq. 7.821.3] and [27, Eq. 9.31.1]. This leads to the following 
closed-form expression for the ergodic capacity of AF MIMO dual-hop systems as the source power p 
grows large for fixed relay power a, 



1. \ 

(66) 



p^OO _ _ 

1=1 k=q-s+l 



4^1 / nsrir 



0, 1, 

k+p + l — q—1, Us + k — q, 0, 

This result shows that if we fix a and take p large, then the ergodic capacity of AF MIMO dual-hop systems 
remains bounded (as a function of a). This confirms the intuitive notion that the capacity is restricted by 
the weakest link in the relay network; in this case, the relay-destination link. 
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Fig. 7. Comparison of bounds, exact analytical, high SNR analytical, and Monte Carlo simulation results 
for ergodic capacity of AF MIMO dual-hop systems with different antenna configurations. Results are 
shown for q//9 = 2. 

V. Tight Bounds on the Ergodic Capacity 

In order to obtain further simplified closed-form results, in this section we derive new upper and lower 
bounds on the ergodic capacity. 



A. Upper Bound 

The following theorem presents a new tight upper bound on the ergodic capacity of AF MIMO dual-hop 
systems. 

Theorem 5: The ergodic capacity of AF MIMO dual-hop systems is upper bounded by 



C{p)<Cu{p) = -log, (/Cdet(3)) 



(67) 



where H is defined in ( [28] ). 



Proof: Application of Jensen's inequality give; 



I 



C{p)^l: log2 E jdet + ^hIlHi 
2 V n.. 



The result now follows by using Theorem |2] 



(68) 



□ 



Note that this inequahty has also been applied in the ergodic capacity analysis of single-user single-hop MIMO systems (see 
eg. [32,38,39]). 



19 

Fig. |7] compares the closed-form upper bound ( [671) with the exact analytical ergodic capacity based on 
( [39l) and ( |40l ). for two different AF MIMO dual-hop system configurations. The results are shown as a 
function of SNR p, with a = 2p. We see that the closed-form upper bound is very tight for all SNRs, for 
both system configurations considered. Moreover, we see that in the low SNR regime (e.g. p « 5 dB), the 
upper bound and exact capacity curves coincide. 

The ensuing corollaries present some example scenarios for which the upper bound (l67l) reduces to 
simplified forms. 

Corollary 6: For the case oo, Cu{p) becomes 

lim Cuip) = Jlog2(/Cdet(Hi)) (69) 
where Hi is a g x g matrix with entries 

{Hi}^ = a'-^'dr-iia) + pa'-^M^) ■ (70) 
Proof: The proof is straightforward and is omitted. □ 
This result shows that in AF MIMO dual-hop systems, when the numbers of antennas at both the relay and 
destination remain fixed, the ergodic capacity remains bounded as the number of source antennas grows 
large. This is in agreement with the results in Section IIV-A. 1 1 

Note that for the scenarios — > cxd and rid —>■ (yo, simplified closed-form results can also be obtained by 
taking the corresponding limits in (l69l) or, alternatively, by using the equivalent single-hop MIMO capacity 
relations in (|42l) and (|44l) . and applying known upper bounds for single-hop MIMO channels in [40]. We 
omit these expressions here for the sake of brevity. 

Corollary 7: Let Ur = 1. Then, Cjj{p) reduces to 

Cij^^Hp) = \^og, (l + /,n,e^i?„,+i (^)) • (71) 
When rid oo, C^''^^{p) becomes 

lim C^^=Hp) = llog^{l + p). (72) 



lim Cjr 



When a ^ oo, C'^'' ^{p) becomes 

lim C^'=i(p) = ilog2(l + p) . (73) 
Proof: See Appendix WD\ □ 
This shows the interesting result that, if a single relay antenna is employed, then when either the number 
of destination antennas or the relay gain a grows large, the ergodic capacity is upper bounded by the 
capacity of an AWGN SISO channel. 



20 



Corollary 8: In the high SNR regime, (i.e. as p ^ cxd) for fixed relay gain a, Cu{p) becomes 

lim Cu{p) = \ log2 f/C det(3)) (74) 
where H is a g x g matrix with entries 



Proof: See Appendix |T 



r(r-l), n<q-ns, 

( \ ^^^^ 

r(T-l) (l + ^ins-q + n)iT-l)), n > q - n^. 

11 □ 



This expression is clearly much simpler than the exact ergodic capacity expression given for this regime 
in (l66l). 



B. Lower Bound 

The following theorem presents a new tight lower bound on the ergodic capacity of AF MIMO dual-hop 
systems. 

Theorem 6: The ergodic capacity of AF MIMO dual-hop systems is lower bounded by 



2 I n.o Is 



^ (n, - s + A:) + /C ^ det (W;. 

fc=l k=q—s+l 



(76) 



where is defined as in (1351) . 

Proof: See Appendix III-FI □ 
In Fig. 121 this closed-form lower bound is compared with the exact ergodic capacity of AF MIMO dual-hop 
systems. Results are shown for different system configurations. The lower bound is clearly seen to be tight 
for the entire range of SNRs. Moreover, in the high SNR regime (e.g. p 15 dB), we see that the lower 
bound and exact capacity curves coincide. 

The ensuing corollaries present some example scenarios for which the lower bound (1761 ) reduces to 
simplified forms. 

Corollary 9: For the case — > cxd, Cl{p) reduces to 

^lirn^ Cl {p) = I log2 [l + pa exp ^ det (W^)^ ^ . (77) 
Proof: See Appendix III-GI □ 
Again, we note that for the scenarios ^ oo and oo, simplified closed-form results can also be 

obtained by taking the corresponding limits in (|69l ) or, alternatively, by using (l42l ) and (l44l ). and applying 
known lower bounds for single -hop MIMO channels in [40]. 
Corollary 10: For the case = 1, Cl{p) reduces to 



CT=' (P) = \ log2 + ;^7^I^ y i^s) + V' (n.) - e(i+^)/° Ei^, 



(78) 
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Fig. 8. Comparison of capacity bounds, high a approximation, and exact analytical results for different 
relay gains. Results are shown for Ur = 1, Us = 2, Ud = 4: and p = lOdB. 

When Us oo, C^"^^ (p) becomes 

^lirn^ ip) = \ log, (^1 + ^ exp {^^ in,) - e^^^^^^ "g' E,^, (i±^) ) ) . (79) 

When Hd oo, C2''^^(p) becomes 

„lim^Cr-(rf = ilog,(l + ;;^»p(«(n,)+*(i±i!))) . (80, 



When a ^ oo, Cl{p) becomes 

lim Cl-=^ [p) = \ log2 + — exp (n,))) . (81) 
Proof: See Appendix III-Hl □ 
As also observed from the upper bound in Corollary [71 this result shows that for a system with a single 
relay antenna, when the relay gain a grows large, the ergodic capacity of an AF MIMO dual-hop channel 
is lower bounded by the capacity of an AWGN SISO channel (with scaled average SNR). 

Fig. E] plots the closed-form upper bound (TtT]) . closed-form lower bound (1781) . and the exact analytical 
ergodic capacity based on ( [39l ) and ( |40l ). for an AF MIMO dual-hop system with = 1. The results 
are presented as a function of the relay gain a. We see that both the upper and lower bounds are quite 
tight for the entire range of a considered. The asymptotic approximations for the upper and lower bounds, 
based on ( [73] ) and ( [STT l respectively, are also shown for further comparison, and are seen to converge for 



22 




Fig. 9. Comparison of capacity bounds, high SNR approximations, and exact analytical results. Results 
are shown for a system configuration (3, 4, 2) and a = 2. 

moderate values of a (e.g. within a 20 dB). 

Corollary 11: In the high SNR regime, (i.e. as p — > cxo) for fixed relay gain a, Cl{p) becomes 

lim Cl{p) = logs I 1 + ^ exp ( - V det (w,) | | , (82) 

\ \ k=q-s+l I I 



where is a g x g matrix with entries 



r - "1 I r (r - 1), n / A; 

{Wfcl =<^ . (83) 

L im,n 1^ r(T-l)[VK-(? + n)+V'(r-l)], n = k 

Proof: See Appendix III-II □ 

As for the high SNR upper bound presented in (1741 ). this closed-form lower bound expression is simpler 

than the exact ergodic capacity expression given for this regime in (l66l l. 

Fig. |9] depicts the closed-form high SNR approximations for the exact ergodic capacity, as well as the 

respective upper and lower bounds, based on (l65l) . (IT?! ), and (l82l ) respectively. For comparison, curves are 

also presented for the upper bound (|67] ). lower bound dTHl ). and the exact analytical ergodic capacity based 

in (l39l ) and (l40l ). Results are shown for an AF MIMO dual-hop system with configuration (3, 4, 2). Clearly, 

the analytical high SNR approximations are seen to be very accurate for even moderate SNR levels (e.g. 

p PS 20 dB). 
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VI. Conclusions 

This paper has presented an analytical characterization of the ergodic capacity of AF MIMO dual-hop 
relay channels under the common assumption that CSI is available at the destination terminal, but not at 
the relay or the source terminal. We derived a new exact expression for the ergodic capacity, as well as 
simplified and insightful closed-form expressions for the high SNR regime. Simplified closed-form upper 
and lower bounds were also presented, which were shown to be tight for all SNRs. The analytical results 
were made possible by first employing random matrix theory techniques to derive new expressions for 
the p.d.f. of an unordered eigenvalue, as well as random determinant results for the equivalent AF MIMO 
dual-hop relay channel, described by a certain product of finite-dimensional complex random matrices. 
The analytical results were validated through comparison with numerical simulations. 
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Appendix I 

Proofs of New Random Matrix Theory Results 

A. Proof of LemmaU} 

To prove this lemma, it is convenient to give a separate treatment for the two cases, q < Ug and q > Ug. 
1) The q < Us Case: For this case, an expression for the p.d.f. has been given previously as 

[25] 

^'"^^^^ ^ Q det (L)--«+i nLi r (n, - z + 1) Ul<j iPj " A) ^^"^^ 
where Z)/ ^ is the k)th cofactor of a q x q matrix with entries 

{i)}^^^=T{ns-q+j)0t~'^' . (85) 

After some basic manipulations, we can express this cofactor as 

6, , = n;..rK-, + i)det(Lr-'^- 

Substituting (l86l ) into ([84b yields the desired result. 
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2) The q > Us Case: For this case, we start by employing a result from [41, Eq. 11] to express the 



joint p.d.f. of the unordered eigenvalues 71, . . . ,7n, of h|lHi, conditioned on L, as follows 



/(7i,--->7njL) 



det(Ai)nr<,(7i-7 



^snr=irK-i+i)m<,(/3i-A 



(87) 



where Ai is the g x g matrix 

1 !3i 



/3 



nq—Us — l TT" 

PI e 



Ai 



The p.d.f. of a single unordered eigenvalue A is found from (1871 ) via 



.q—ris — l 



e "1 



(88) 



/a|l(A) 



CO roo 



/ (71, . . . , 7„J L) (i7i • • • d7n,-i 



1 



ns nr=i r{ns-i + l) UUj Wj - P^J JO 



7n.=A 
00 />oo 



det(Ai)det(7f Md7i • • • d7„^_i 



7n.=A 

(89) 



where we have used Y[7<j (^i ~ 7«) ~ '^^^ l^^^j' evaluate the Ug — 1 integrals, we expand det (Ai) 
along its last column and det (^7i along its last row, and then integrate term-by-term by virtue of [33, 
Lemma 2]. This yields 

q <? 



/ail (A) 



1=1 k=q—ns+l 



ri.nr=irK-i+i)n?<,(/?i-/3.) 

where Di^k is the (/, k)th cofactor of a g x g matrix H = [AC], with entries 



(90) 



{A}m,n = /5m^ m=l,...,q, n=l,...,g-n. 



(91) 



and 



{cUn = rH/3; 



.q—ris+n—l 



m = 1, . . . ,q, n = 1, 



(92) 



Then, it can be shown that 

q Q 



1=1 k=q—ns+l 



k=q—ns+l 



(93) 
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where Dy^ is a q x q matrix with entries 



rn = l 

r {n - q + Us - 1) i3^, m = 1 
/3«r"'"'e-V/3™A"-«+"=-i, m=l,...,q, 



E det (Dfc) 



Hence, we can rewrite (1901 ) as follows 



,...,q, n = l,...,q-ns 
,...,q, n = q - Us + 1, . . . ,q, n^k 



n = k 



/ail (A) 



n.nr=ir(n.-i+i)n?<,(/3.-A 

After some basic manipulations, ( [95] ) can be further simplified as 

1 



/a|l (A) 



1 ^n,-q+k-l 

E fT:^ det (D,) 



where is a g x (7 matrix with entries 



I. ""J m,r. 



Pm 1 n k, 

e-V/3™/3^-"=+\ n = k. 

Finally, we apply Laplace's expansion to ( l96l ) to yield the desired result. 



(94) 



(95) 



(96) 



(97) 



B. Proof of Lemma^ 

The joint p.d.f. of Wi = diag {ai, . . . , aq} is given by [42-44] 



/Wi ("1, 



/Ce'£"narn(a, 



i=l 



i<j 



Recalling that 



ai 



Pi 



(98) 



(99) 



1 — aP, 

we derive the joint p.d.f. of W2 = diag{/3i, . . . from ( [98] ) by applying a vector transformation [45] 

Pi Pq 



/ws {Pir-- ,Pq) = /Wi 



1 — a/5i ' ' 1 — a/3q 



\3{{ai,...,ag) ^ (100) 



where 



J ((ai, ...,aq)^ (Pi, . . .,Pq)) = det 



9a 1 

a/3i 



S/3i 



8a 1 



(101) 
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From (|99l ). we have 

dai 1 

" (l-aft)2' 

therefore the Jacobian transformation in (IIOII ) is evaluated as 

J ((«!,..., ag) ^ (/?!,..., =TT- — 2- 

f=i (1 - aA) 

Substituting ^ and (fT03] ) into (fTOOl ) yields 

|( 1 - a A )P~ +2 J- i V 1 - 1 - a A 

Finally, simplifying using 

I/O r, \ 2 Q 



(102) 



(103) 



(104) 



TT r - V = TT { V 



(105) 



yields the joint p.d.f. of L. 

We now derive the p.d.f. of an unordered eigenvalue /? of the diagonal matrix L. According to [31, Eq 
42], the unordered eigenvalue p.d.f. of H2H2 is given by 

q~\ i 2j 

q 



/ (^) = - E E E ^ ^' A^^'^+'e-^ . (106) 



i=0 j=0 1=0 

Recalling that /? = A/ (1 + aX), the result follows after applying a simple transformation. 
C. Proof of Theorem [7] 

We start by re-expressing the conditional unordered eigenvalue p.d.f. in Lemma\l}as follows 

^ Jn^TM '''' H • 

where Dfc is a g x g matrix with entries 

=\ ^ (108) 

\ e-V/3™/3^-"=-i, n = k. 

Now, utilizing Lemma El we can evaluate the unconditional p.d.f. as follows 

/A(A)=ii;L[/A|L(A)] 



ir 1 \ns~g+fc— 1 
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where 

1 



r -/ y op-q -jr^ 

/ det(D,)n(/3,-A)n ri aB)pW 



• • • 



= det(Yfc), (110) 
where is a. q x q matrix with entries 

{fl/a ».p-g+»n+7i-2 S — 

Let t = x/(l-ax). Utilizing [27, Eq. 3.383.5] and [27, Eq. 3.471.9], the integrals in (fTTTT) can be 
evaluated, respectively, asc 

^ —e-T^dx= / (1 + at)^^-™-" e-*di 

(1 - axY^^ Jo 

= a9-P-'^-"+ir {p - q + m + n - 1)U {p - q + m + n - l,p + q,l/a) 

(112) 



and 



l/a p— ria+m— 2 



(1 - ax)P+^ 
= e 

g+ria— m 



Jo 



j=0 

where U (•, •, •) is the confluent hypergeometric function of the second kind [27, Eq. 9.211.4]. 
Combining (|109l )- (ll 131 ) and then applying Laplace's expansion yields the desired result. 

D. Proof of Lemma \3\ 

We will prove the lemma by giving a separate treatment for the two cases, q < rig and q > Hs 



Note that, by using the Binomial expansion, ill2t can be alternatively expressed as 

/ t''-''+"'+"--\l + atf e-'dt= V a'r{p~q + m + n + i-l) 

J j-n 
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1} q < Us Case: In this case, we start by writing 



^ <; det ( + — h|lHi 



L det ( Ig + — LHiH| 



^ n 1+ 



.1=1 



pa 



-7! 



(114) 



where 71, . . . , 7^ are the ordered eigenvalues of LHiH|. Conditioned on L, the joint p.d.f. of 71, . . . , 7^ 
is given in [46]. Using this result, we can express (II 141 ) as follows 



^ <! det ( + — h|lHi 



L det (e-^^V/^O nLi 1 + ^70 /5r"=-Sr="^ det(7r^)ci7i • • • ^7, 



(115) 



nLirK-i + i)n?<,(/?,-A) 

where the integrals are taken over the region Pord = {c« > 7i > • • • 7ij > 0}. Applying [46, Corollary 2], 
(II 151 ) can be evaluated in closed-form as follows 

nLi/3r"="'det(Hi) 



s <; det { In, + — hIlhi 



nLir(n.-i + i)m<,-(/5.--A)' 



(116) 



where Hi is a g x g matrix with entries 



{Hi}^^„ = (^r (n, - g + n) + ^/3„r (n, - g + n + 1)^ (1 17) 



Extracting common factors from the determinant in (II 161 ) and simplifying yields the desired result. 
2j q > rig Case: In this case, we use the joint eigenvalue p.d.f. (l87l) to obtain 



^ <; det ( I„, + ^h|lHi 



Tig 



L =E n 1+ 



.1=1 



pa 



rig 



nr=i ( 1 + f 7.) det (Ai) det(7r^)d7i • • • ^7 



(118) 



where 71, ... , 7„^ are the ordered eigenvalues of h|lHi, Ai is defined in ( |88l ). and the integration region 
is Pord = {00 > 71 > • • • 7n^ > 0}. Applying [33, Lemma 2], dl 181 ) can evaluated in closed-form as 
follows 



<{ det { + — hJlHi 



Ug 



det (H2) 



(119) 



where H2 = [Ai Ci] is a g x g matrix with entries 



{Ai}™„ = /?;^ , n = l,...,g-re^ 



(120) 
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and 

{Ci}^,„ = C+'^-"^-'(r(n) + (/>aM)/3„r(n + l)), n=l,...,n,. (121) 
Extracting common factors from det (H2) and simplifying yields the desired result. 

E. Proof of Lemma^ 

To prove this lemma, it is convenient give a separate treatment for the two cases, q < Ug and q > Ug. 
1) q < Us Case: Now we need to calculate the expectation E jlndet (^LHiHJ)}. The moment 
generating function (m.g.f.) of In det (LHiH| j, conditioned on L, is given by 



. (122) 



Mi{t\L) = E^det [LUiH[ 
Utilizing the joint p.d.f. of the eigenvalues 71, ... ,7^ of LHiH|, presented in [25,46], we get 

where the integrals are taken over the region J^ord = {00 >7i > •••7^ >0}. Applying [46, Corollary 2], 
(11231 ) can be further simplified as follows 

\\i=i T{ng-i + l) ni<j - A) 

where 33 is a g x g matrix with entries 

poo 

{33}^ n = / e-^/^"y--«+*+"-idy = Pt^^'^T (n, - g + t + n) . (125) 

From M.I {t |L), we get 

E^jlndet (^LHiHJ) L}=^A^i(t|L) 



E det (Sfc) 

k=l 



nLirK-i+i)nr<,(/3i-A 

where is a g x g matrix whose entries are 



(126) 



(hs- q + n) {us - q + n) + ln/?m], n = k. 

where ^(•) is the digamma function. Now, det (Xlfc) can be further simplified as 

g 

det (Sfc) = det J{V{ns-q + k) (128) 

k=l 
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where is a g x g matrix with entries 



{^4 



gn-1 



n ^ k, 



(129) 



^['ip{ns-q + n)+ In Prn], n = k. 



By using the multi-linear property of determinants, along with some basic manipulations, we can write 



det ( Sfc ) =ijj{ns-q + k) det ( ) + det (Y 



(130) 



Substituting (11281 ) and ( 11301 ) into ( 11261 ) and simplifying yields the desired result. 

2j q > rig Case: We now evaluate the m.g.f. of In det (^hIlHi^ conditioned on L, which is given by 

(131) 



M2{t\'L) = E\ det H{LH 



Utilizing dSTl ). (11311 ) can be expressed as 

1 



M2{t\L) 



nr=i r (n, - z + 1) ut<j - ft) t\ 

where Pord = {oo > 71 > • • • jn, > 0}. Applying [33, Lemma 2] yields 

det (H4) 



/ rr7*det(A2)det(7/-i)d7i,...,d7„^, (132) 



M2{t\'L) 



nr=ir(n.-i+i)n?<,(/3,-ft)' 

where 34 = [A2 C2] is a q x q matrix with entries 



(133) 



{^2}^,n = Pm \ n= l,...,q 



(134) 



and 



{C2}„,„ = r(t + n)/?^- 



iq—ris+t+n—l 



n= 1, . . . ,n<j 



(135) 



From the m.g.f. (11331) . we can then obtain 

E^jlndet (^h|lHi) L}=^A^2(i|L) 



E det (ftk) 

k=q—ns+l 



nr=irK-i+i)n?<,(/?i-ft 



(136) 
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where fi^ is a g x g matrix with entries 

r'm ) 



L J m,n 



T{ns-q + n)P:^-\ 



n^k, n = I, . . . ,q - Us, 
n^k, n = q - Us + 1, . . . ,qs, 



(31^ T {us- q + n) [V" [ug - q + n) + In/?™], n = k. 



(137) 



By using the multi-Unear property of determinants, along with some basic manipulations, we can obtain 
the desired result. 

3) q = s Case: In this case, starting with (|3TI) . we can write the determinant summation over k as 
follows 



^det (Yfe) = ^^sgn(a) 

k=l {a} 



k=l 



.i=l 



ln/3„(fc) 



^det (Yfc) = ^sgn(a) 

{a} 



k=l 



•a(k) 



.4=1 J fc=l 



lndet(diag{A}ti)n-.(^^-/^' 



lndet(L)rr' .iPj-Pi) ■ 



Substituting (11391 ) into (|3TI ) yields the final result. 



(138) 



where the second summation is over all permutations a = {a (1) , . . . , a (g)} of the set {1, . . . ,q}, with 
sgn(a) denoting the sign of the permutation. We can further write 



(139) 



F. Proof of Theorem \3\ 

We start with Lemma E] and remove the conditioning on L by using Lemma \2\ as follows 

s 

E {In det (*)} = ^1!^ {hs - s + k) 



k=l 



0</3i<-<A 



det (/^r'j n 5 (A) Yl det (Y,.)(i/?i ■■■d(3g, 

'J^^/'^ i=l fc=g-n,+l 



where 



yP-q^-u/(l-au) 



(1 - au)P+^ 

Using [33, Lemma 2], these integrals can be simplified to give 

s q ^ 

E{lndet{^)} = Y^ins-s + k) + IC ^ det (w^ 

k=l k=q—ns+l 



(140) 



(141) 



(142) 
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where is a q x q matrix with entries 

{w,} = 



fo^'' ""i^anK+^ e 1—' Inudu, n = k. 



(143) 



For the case n k, a closed-form expression is given in (II 12I ). For the case n = A;, we utilize [27, Eq. 
4.358.5] and [31, Eq. 47], to obtain 



l/a ^p—q+m+n—2 



e 1-"" Inudu 



(1 - auf^'i 

roo 

= / tP-a+"^+"-2 (1 + at)29-"^-« e-t [i^ t - In (1 + at)] dt 
Jo 



2q—m—n 



j=0 
2g— m~n 



i=0 



^p+q-i-2^-t [int _ In (1 + at)] 
T {p + q - i - I) 

p+q—i-2 



1=0 



(144) 



Substituting (11121 ) and (11441 ) into (1143b and (11421) yields (l33l) . 

When g = s, we start with (l33l) and remove the conditioning on L as follows 



E {lndet{^)} = Y^{ns-q + k) + q / f {p)ln$dp 



(145) 



where / denotes the unordered eigenvalue p.d.f. of L (i.e. p.d.f. of a randomly-selected /9 G {/3i, • • • ,(3q}). 
Substituting this p.d.f. from ([TS] ) and integrating using (11441 ). we obtain the desired result. 

Appendix II 
Ergodic Capacity Proofs 

A. Proof of Eq. dH]) 

When iir oo, the ergodic capacity expression ([TT]) can be expressed as follows 



1 



lim C{p) = -^<^log2det + 



2 t °^ V (1 + p) ' 

where Li = diag {A?/ (n^. (l + aA^)) }. Noting that q = Ud, by the Law of Large Numbers we have 



(146) 



lim 



H2Ht 



Ur^OO fir 



(147) 



which implies that 



lim 



rir-^oo fir 



1 , i = l,...,nd ■ 



(148) 
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Recalling application of ( 11481 ) in ( 11461 ) yields 



lim C ip) = -E log2 det I,^ + ^ H^H , (149) 

ur-^oa 2 [ \ ns(l + p + a) J] 

where H is an x ?is i.i.d. Rayleigh fading MIMO channel matrix. Applying the identity ^ to ( |149b 
yields the desired result. 



B. Proof of Eq. (EJ 

Using ([6]), the ergodic capacity expression ([TT]) can be alternatively written as 

C (p) = ^E |log2 det (^Ig + ^HiHIl^ I . (150) 

By the Law of Large Numbers we have 

HiH| 

lim I„ (151) 



and hence (11501 ) reduces to 



lim C (p) = -E {Iog2 det (Ig + paL)} . (152) 



Substituting ([121 ) into ( 11521 ). after some simple manipulations we easily obtain 

Jlrn^ C (p) = ]^E |log2 det (l^ + (p + 1) aYlln^ } {l°S2 ^et (ig + aH^H2) } . (153) 
Substituting ([5]) into (11531 ) and applying the identity ^ yields the desired result. 

C. Proof of Theorem [?] 
We will consider the following cases separately; namely, q < Ug and q > Ug. 

1) q < Us Case: We start by applying the identity (O to obtain the ergodic capacity, in the high SNR 
regime, as follows 

1 



C{p)\ 



ci,p— »oo,a/p=/3 2 



qlog2 p - (?log2 ( ^ ) + E {log2 det (lHiHO | 



(154) 



The high SNR slope can be calculated as 

5oo = f bit/s/Hz(3dB) • (155) 
Applying ( |49l ). the high SNR power offset is given by 

= 1 (;^) - {^°^^ (^^^^0 } • ^'''^ 
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Invoking Theorem\3\ and simplifying yields the high SNR power offset for case q < Ug. 

The proof of ( [53] ) follows along similar lines to that used above, but in this case invoking Theorem \3\ 
in place of Theorem^ 

2) q > rig Case: In the high SNR regime, the ergodic capacity can be approximated as 



Cip)\ 



_ 1 



ris log2 (p) - ns log2 



+ ^{log2det (hIlHi)} 



In this case, the high SNR slope is 



71 

Soo = ^ bits/s/Hz (3dB) 



and the high SNR power offset can be obtained as 
The result follows by applying Theorem |21 



(157) 



(158) 



(159) 



D. Proof of Corollary [7| 

Substituting = 1 into (l67l l yields 



1 



V nd,nd + 1: 



1 + P 



a 



+ pndU [ nd + l,nd + I 



a 



Using the following properties of the confluent hypergeometric function of the second kind [27] : 



(160) 



U{a,a,z) = e'z'-''Ea{z) 



(161) 



and 



U{a,a + l,z) = z'"", 



(162) 



we get the final expression for C^'-^^p) in ([IB- Note that C^=^(p) can be lower and upper bounded as 



n,,=l / 



(163) 



with 



1 



+ nd+l 



(164) 



and 



CuyHp) = i^^ogJl + pnd- 



(165) 
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where we have used the inequaUty [36, Eq. 5.1.19]. Taking — > oo, we see that both (11641 ) and (11651 ) 
converge to the same limit in dTll) . Taking a ^ oo and ultilizing [36, Eq. 5.1.23], we obtain (173] ). 

E. Proof of Corollary |S] 

Note that when p ^ oo, then a — > 0. Therefore, we apply the following asymptotic first-order expansion 
for the confluent hypergeometric function [36] 

C/(c,6,z) = z-^ + o(l), z^oo (166) 

to yield the desired result. 

F. Proof of Theorem |6| 

We will use the lower bound derived in [40, Theorem 1] and consider the following cases separately; 
namely, q < Ug and q > Us. 

1) q < Ug Case: Applying the ^ and [40, Theorem 1] to (ITTI ). we lower bound the ergodic capacity, 
conditioned on L, as follows 

C(p) > glog2 (^1 + ;^exp Q^jlndet (^LHiH|)}^^ . (167) 

Now, using Theorem \3\ yields the desired result. 

2) q > rig Case: In this case, the lower bound can be written as 

C(p)>n,log2(^l + -^exp(^-^-E|lndet(^HiLHl)}^^ . (168) 
Again, we use Theorem \3] to obtain the desired result. 

G. Proof of Corollary |9] 

When Ug —>■ oo, ip [ug — q + k) can be approximated as [36, Eq. 6.3.18] 

ilj{ns-q + /c)|„^^oo ~ In (n^ - ^ + k) 

K.\ang . (169) 

Substituting (11691 ) into (1761 ) yields the desired result. 

H. Proof of Corollary \M 

Taking ^ oo and using [36, Eq. 6.3.18], we get ( [79l) . 
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For the case — > oo, we first apply [36, Eq. 5.1.19] and [27, Eq. 8.365.3] to obtain the following 
approximation 

Furthermore, substituting (11701 ) into ^ and using [27, Eq. 8.365.5] and [36, Eq. 6.3.18] yields (l80l ). 

Now consider the case a ^ cxd. Utilizing the recurrence relation for the exponential integral [36, Eq. 
5.1.14], the summation in (iTSl) can be alternatively written as 




(171) 



where 7 = 0.577215 ... is the Euler's constant. Note that, in deriving ( 11711 ). we have applied the definition 
of the digamma function [27, Eq. 8.365.4]. Using the series expansion given in [36, Eq. 5.1.11], when 

a ^ 00, we get 



and therefore 




(172) 



(173) 



Applying (I171I )- (I173I ) in (iTSl ) yields the desired result. 

/. Proof of Corollary |77] 

Using the following approximation [36] 

^,(z)^ie-^(l + o(i)) |z|^oo, (174) 

qm+n{a) can be approximated as 

W+n(a)|,_.^ « r (r - 1) V (r - 1) , (175) 

which leads to the final result. 
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